Conversed spin Hall conductance in two dimensional electron gas in a perpendicular 

magnetic field 
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Using the microscopic theory of the conserved spin current [Phys. Rev. Lett. 96, 076604 (2006)], 
we investigate the spin Hall effect in the two dimensional electron gas system with a perpendicular 
magnetic field. The spin Hall conductance <r^„ as a response to the electric field consists of two 
parts, i.e., the conventional part <r*° and the spin torque dipole correction a s J v . It is shown that the 
spin-orbit coupling competes with Zeeman splitting by introducing additional degeneracies between 
different Landau levels at certain values of magnetic field. These degeneracies, if occurring at 
the Fermi level, turn to give rise to resonances in both <t^° and afZ in spin Hall conductance. 
Remarkably, both of these two components have the same sign in the wide range of variation in the 
magnetic field, which result in an overall enhancement of the total spin Hall current. In particular, 
the magnitude of a s J v is much larger than that of crf° around the resonance. 



PACS numbers: 72.25.Dc, 72.10.Bg, 73.23.-b 
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Spintronics, which combines the basic quantum me- 
chanics of coherent spin dynamics and technologi- 
cal applications in information processing and storage 
devicesi^i 3 -, has grown up to become a very active field in 
condensed matter physics. One central issue in spintron- 
ics is on how to generate and manipulate spin current as 
well as to exploit its various effects. In the ideal situa- 
tion where spin is a good quantum number, spin current 
is simply defined as the difference between the currents 
of electron carried by spin-up and spin-down states. This 
concept of the spin current has served well in early stud- 
ies of spin-dependent transport in metals and although 
the ubiquitous presence of spin-orbit coupling inevitably 
makes the spin non-conserved, this treatment on spin 
current sustains to be reasonable if the net effect of spin- 
orbit coupling is only considered as one source to spin re- 
laxation. In recent years, however, it has been found that 
the extrinsic or intrinsic spin-orbit coupling can provide 
a route to generate transverse spin current in ferromag- 
netic metals^ or semiconductor paramagnets^ 7 -^ by the 
driving of an electric field. The fundamental question 
of how to define the spin current properly in the gen- 
eral situation then needs to be answere d 9 ' 10 ! 11 ' 12 ! 13 ' 14 . In 
most of previous studies of bulk spin transport, it has 
been conventional to define the spin current simply as a 
combined thermodynamic and quantum-mechanical av- 
erage over the symmetric product of spin and velocity 
operators. Unfortunately, the recent spin- accumulation 
experiment a 15 ' 16 ' 17 cannot directly verify it since there is 
no deterministic relation between this spin current and 
the boundary spin accumulation. 

Recently, Shi et al3 have established an alternative 
definition of spin current, which is given by the time 
derivative of the spin displacement (product of spin 
and position observable). This differs from the con- 
ventional definition by inclusion of one additional spin 
torque dipole term. As a result, the new spin trans- 
port coefficients a* have been shown to consist of two 
parts, i.e., the conventional part afj^ and the spin torque 



dipole correction part a s ^ v . As one key consequence, the 
Onsager relation between <r* and other kinds of force- 
driven transport coefficients has been shown 9 -^. A gen- 
eral Kubo formula for the spin transport coefficients a s ^ v 
in terms of single-particle Bloch states has been giver* 9 ^ 
and applied to study the conserved spin Hall conductiv- 
ity (SHC) in the two dimensional hole gas (2DHG). Also, 
the SHC based on the new spin current has been recently 
calculated 2 ^ for two dimensional electron gas (2DEG) 
with /c-linear or fc-cubic spin-orbit coupling and found 
to depend explicitly on the scattering potentials. 

In this paper we study the conserved spin Hall effect of 
2DEG with Rashba spin-orbit coupling in a perpendicu- 
lar magnetic field. This system with the same setup has 
been studied recently by Shen et ali^i within the conven- 
tional spin current framework. They have found that the 
conventional SHC can be made resonant or even diver- 
gent by turning the sample parameters and/or magnetic 
field B. The behavior of the spin torque dipole contribu- 
tion, and subsequently the total conserved SHC remains 
yet to be exploited. As we will show in this paper, the 
torque term in the SHC can also be made resonant or di- 
vergent by tuning the magnetic field or the Rashba spin 
splitting. Moreover, we find that the resonant ampli- 
tude of torque term <j s J v is even more prominent than 
that of the conventional term o s ® v . The oscillations and 
resonances in both of these two terms stem from energy 
crossing of different Landau levels near the Fermi level 
due to the competition between Zeeman energy splitting 
and spin-orbit coupling. Another fact we find is that dif- 
ferent from the case without magnetic field, the spin Hall 
currents contributed from the conventional part and the 
spin torque dipole part flow in the same direction, which 
means an overall enhancement of the total conserved spin 
Hall current compared to the conventional spin Hall cur- 
rent. 

We consider a 2DEG with the Rashba coupling in the 
x-y plane of an area L x L subject to a perpendicular 
magnetic filed B = —Bz. The electrons are confined 
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FIG. 1: (Color online). Illustration of the two-dimensional 
electron system in the combined magnetic and electric fields. 



between —L/2 and L/2 in the y direction by an infinite 
potential wall, and its wave function is periodic along 
the x direction. We choose the Landau gauge A = xBy. 
The Hamiltonian for a single electron of spin- 1/2 with a 
Rashba coupling is given by 

II 2 A 1 

H ° = 2m. + h^ x<Jy ~ n * <7 ' t ) ~ 2 9s ^ bB(Jzi D 

where the confining potential is implied, m, (— e), and 
g s are the electron's effective mass, charge and effective 
magnetic factor respectively, fis is the Bohr magneton, 
II = p + eA/c is the kinetic operator, A is the Rashba 
coupling, and a a are the Pauli matrices. In the presence 
of the Rashba spin-orbit coupling, the electron momen- 
tum p x = hk along the x-direction remains to be a good 
quantum number, thus in the Hilbert subspace of given 
k, the Hamiltonian can be written as 



Ho(k) = fko 
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where uj = eB/mc is the cyclotron frequency, r\ — 
Xmlb/h 2 is the effective Rashba coupling, and g — 
g s m/2m e , with m e the mass of a free electron and 
h = \fhc/eB the magnetic length. a& = [y + (fifc + 

ipy)c/eB]/V2lb, so that [<Zfc,af.] = 1, a± = (a x ± i<J v )/2 
are the Pauli matrices. The energy spectrum of Ho(k) is 
given by 



e ns = hj \ n + - v(l 



g) 2 + 8nrj- 



(3) 



with s = ±1, for n > 1; and s — 1 for n = 0. The cor- 
responding two-component eigenstates for e ns are given 
by 



\nks) = 



cos ns \n)k 
isin8 ns \n — l)k 



(4) 



where \n)k is the eigenstate of the nth Landau level in the 
absence of the Rashba interaction. For n = 0, 6*oi = 0, 
otherwise for n > 1, tan6>„ s = — u n + + u 2 , with 
u n = (1 — g)/\/8nr]. The eigenstate |n, k,s) has a de- 
generacy — L 2 eB/hc, corresponding to quantum 
values of k. 

Now let us consider a uniform electric field E applied 
in the y-direction. The total Hamiltonian in this case is 
H = Ho+eEy in which the term eEy is usually treated as 



a small perturbation. Within the conserved spin current 
formalism, the spin current is defined as a time-derivative 
of the spin displacement operator, i.e., J s = d(is z )/dt, 
where s z is the spin operator for a particular component 
(z here, to be specific), and f is the electron position op- 
erator. Compared to the conventional spin current oper- 
ator, there has an extra term, f (ds z /dt) in J~ s , which ac- 
counts the contribution from the spin torque dipole. As a 
consequence, the conserved SHC crj as a linear response 
to the external electric field defined by S s , x = VxyE, 
includes two components, 

a xy — a xy ' a xy \° ) 

The first term in Eq. ([5]) is the usual conventional SHC, 
which is ready to be rewritten in the Landau spinor space 
(consisting of states \nks)) as 

eh \ - 
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[f(ens)-f(en's>)]X (6) 



lm(nks\^{v x , s z }\n' ks') (n' ks' \ v v \nks) 

(tnks ~ tn'ks') 2 + S 2 

where the velocities are given by 
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and f(e ns ) is the equilibrium fermi function. The limit 
of 8 — > in Eq. © should be taken at the last step 
of calculation, and as a result, there is no intra-band 
({n, s} = {n',s'}) contribution. The Kubo formula (|6|) 
for the conventional SHC has been used by most of pre- 
vious investigations. The second component a x T y in the 
conserved SHC (5) comes from the contribution of the 
spin torque dipole terrn^. In the present context, its ex- 
pression reads 



st efl r 1 

y V q^O q 



^2 if( e nks) ~ f(e n 'k+qs')] X 



(8) 



Rc[(nfcs|f (fc, q)\n'k + qs')(n'k + qs'\v y (k, q)\nks)\ 

(fnfcs — Cn'fc+gs') 2 + ^ 

where T s (k,q) = \ [T s (k) +T s (k + q)], v y (k,q) = \ [v y {k) + 
v y (k + q)], with t s = j^[s z ,Hq\. Note that to properly 
calculate a a J in practice, all the terms in Eq. ([5]) with 
the subscript k + q should be expanded at k to first order 
in q. 

By substitution of the expressions for single-particle 
velocities v x , y and eigenstates \nks), and after a tedious 
but straightforward derivation, we obtain the conven- 
tional part of conserved SHC as follows 
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where 



A(n, s, a') = (n + 1) cos 2 9 ns cos 2 9 n+ls , (10) 



nsin 2 9 ns sin 2 Q n+ i s > 



+ -j=(Vn + lcos 2 9 ns sin26»„ + i s / 

- y/nsin 2 9 n+ \ S ' sin2# ras ). 

The other component contributed from the spin torque 
dipole term turns out to be given by 



2e?7 f( e ns) - f(e n+ w) 



ST 

u *v = — Is 

2er) 



n,s,s' 



B{n,s,s') (11) 



where 



B(n,a,a') = [ncos9 ns sm9 n+ i^ 

- (n + 1) sin 6> n;5 cos 0„+i jS ' 



(12) 



x [ ^/nfn + 1) cos ns cos 9 n +\s 
+ nsin 9 ns sin0„ + iy 

+ \ph\r\ cos 9 ns sin n+ \ s i] 



and 



C(n, s, s ) = sin (9., 



Ons') x [n 2 cos9 ns cos9 ns > (13) 
+ (n - 1) \/n(n - 1) sin 9 ns sin ns / 



-\-n\l1nr) cos s sin # ns / ] . 

A similar expression for the conventional component 
a^y has also been obtained by Shen et ali^L, while the 
spin torque dipole term u 8 J y is derived for the first time 
in this paper. The difference between <j s ® y and g s J v is 
obvious. Note that the first line in Eq. (fTTf is obtained 
by expanding \n'k + qs'} in Eq. © at k to the first order 
in q while remaining other quantities to be their values at 
q = 0. Whereas the second line in Eq. (|11|) is obtained 
by a linear expansion of (n'k + qs'\ with respect to q. 
The other terms occurring in Eq. ([8]) turn out to take no 
contributions to o s J y for the present model. 

From Eqs. and (fTTj) , specially at zero temperature, 
one can see that if \n, s) and \n+ 1, a') are both occupied 
or empty, then the contribution from these two states to 
the SHC is zero. Thus only the states near the Fermi 
level are important to the SHC; this happens in the situ- 
ation that the energy-lower state \n, +1) is fully occupied 
while the upper state |n + l, — 1) is partially occupied or 
fully empty. Due to the denominator in Eqs. (|9|) and 
(fTTj) . one can find that there will occur resonances at the 
degeneracies between the Landau levels e Ul i and e n +i,-i- 
To show this, we first plot in Fig. 2 the Landau lev- 
els as a function of the effective Rashba coefficient 77 for 
experimentally accessible Ino.53Gao.47As/Ino.52Gao.4sAs 
system 22 . Here the effective magnetic factor is taken to 
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FIG. 2: Electron Landau levels as a function of effective 
Rashba coupling rj = Xmlb/h 2 for g=0.1. 
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FIG. 3: (Color online). Conserved SHC and its two compo- 
nents versus 1/B at T = 0. The parameters are A = 9 meV 
nm, n e = 1.9 x 10 _2 /nm 2 , g s = 4, and m = 0.05m e , taken for 
the inversion heterostructure Ino.53Gao.47As/Ino.52Alo.4sAs. 



be g = g s m/2m e = 0.1. One can see that with increasing 
the amplitude of the effective Rashba coefficient rj, the 
pair of states |n, a = 1) and \n + l,s' = —1) approach 
to cross, implying resonances in conserved SHC at these 
level crossings. From expressions for Landau levels e n ,s, 
one can see that the resonant condition for conserved 
SHC is given by 



y/(l - g) 2 + 8m7 2 + y/(l - gf + 8(n + l)r, 2 = 2, (14) 

where 2n < v < 2n+ 1, n = 0, 1, 2, • • • , and v is the filling 
factor given by the ratio of the total electron number 
= Snks /( f nks) to the Landau level degeneracy factor 
Nt/,, i.e., v = Ne/Nt/,. Eq.fJJJ ensures that the state \n, 1) 
is fully occupied, while the state |n + 1, —1) is not fully 
filled. Based on Eqs. (|9))- (fT3"|) we have systematically 
calculated the SHC in a wide range of system parameters. 
Figure 3 plots o s xy and its two components as a function of 
inverse of magnetic field (or the effective Rashba coupling 
coefficient r\). One can see that there is a pronounced 
resonance at filling v = 12.6. At this filling the 13th 
Landau level (state \n = 7,-1)) is partially occupied, 
while the 12th Landau level (state \n = 6,1)) is fully 
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FIG. 4: (Color online). Conserved SHC and its two compo- 
nents versus Rashba coupling A at T = 0. The parameters 
are n e = 1.9 x 10 _2 /nm 2 , g s = 4, and m = 0.05m e . 

occupied. At the same time these two levels cross at 
v = 12.6, as shown in Fig. 2, therefore resulting in a giant 
resonance revealed in Fig. 3. Note that the resonant 
conditions are the same for the conventional term a s ® y 
and the spin torque dipole term cr%L, so one can find 
in Fig. 3 that the resonant peaks in a s xy and a^ y have 
the same location (at B = 6 T). Furthermore, it reveals 
in Fig. 3 that these two conductance components have 
the same sign in a wide range of magnetic field, and the 
resonant amplitude of a s J y is even more larger than that of 
a^y. As a result the total conserved SHC is prominently 
enhanced by the inclusion of the spin torque dipole term 
in. This result is different from that in the absence of 
the magnetic field, wherein it has been verified that the 
two components cr*° and a%L always compete each other, 
giving the opposite contributions to the total spin Hall 
current^. Note that the small side peaks occurred in Fig. 



3 reflect the variations in filling factor v as a function of 
the magnetic field. 

We also calculate the conserved SHC by only varying 
the Rashba coefficient A while the magnetic field remains 
unchanged. The result is shown in Fig. 4 for B = 6 
T. The other parameters are set to be the same as used 
in Fig. 3. One can see that at A = 9 meV nm, there 
occurs a resonance in both a^y and a 8 J y with different 
peak amplitudes and line widths, the resonance behavior 
is even more prominent for u s J v . As a consequence, the 
resonant features of total conserved SHC is dominated 
by its spin torque dipole term. The resonance shown in 
Fig. 4 corresponds to the same case as Fig. 3 does, i.e., 
the 12th Landau level is fully occupied while the 13th 
Landau level is not fully filled. The difference is that 
in Fig. 4 the filling factor v remains unchanged when 
varying the Rashba coefficient. As a consequence, unlike 
what is shown Fig. 3, it reveals in Fig. 4 that there are 
no side peaks occurred. 

In summary, we have studied the spin Hall effect in 
2DEG system in a perpendicular magnetic field by us- 
ing the conserved definition of spin current, which in- 
cludes both the conventional part <r*° and the spin torque 
dipole correction term a 8 J v . We have shown that the con- 
served SHC can be featured by giant resonances by tun- 
ing the amplitude of system parameters. Furthermore, it 
has been found that compared to the previous result of 
cr^y— , the resonance features, including the height and 
line- width of the resonant peak, are even more prominent 
for o s J y . It is expected that the present results could have 
helpful implications on other aspects involving spin trans- 
port and spin accumulation in the presence of magnetic 
field. 

This work was supported by CNSF under Grant No. 
10544004 and 10604010. 
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